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A POSTERIORI ERROR ESTIMATES FOR DISCONTINUOUS 
GALERKIN METHOD TO THE ELASTICITY PROBLEM 

LUONG THI HONG CAM AND DAVEAU CHRISTIAN 


Abstract. This work concerns with the discontinuous Galerkin (DG) method 
for the time-dependent linear elasticity problem. We derive the a posteriori 
error bounds for semi-discrete and fully discrete problems, by making use of 
the stationary elasticity reconstruction technique which allows to estimate 
the error for time-dependent problem through the error estimation of the 
associated stationary elasticity problem. To this end, to derive the error 
bound for the stationary problem, we present two methods to obtain two 
different a posteriori error bounds, by L 2 duality technique and via energy 
norm. For fully discrete scheme, we make use of the backward-Euler scheme 
and an appropriate space-time reconstruction. The technique here can be 
applicable for a variety of DG methods as well. 


Introduction 

A posteriori error estimation and adaptivity recently have become success¬ 
ful tools for efficient numerical computations. Many adaptive methods have 
already been widely considered for solving elliptic, parabolic and first order hy¬ 
perbolic problems. For example, in jT] they consider the a posteriori error esti¬ 
mates for Maxwell equation, [2] considers a posteriori for time dependent Stokes 
equation, [3] investigates the a posteriori error estimate for parabolic problem. 
Howerver, there are few results about a posteriori error analysis for the second 
order evolution problems; we mention, in particular, [3j and |5j derive rigor¬ 
ous a posteriori bounds for conforming finite element methods in case of fully 
discretization for wave equation. 

An important class of finite element methods goes under the general term dis¬ 
continuous Galerkin (DG) methods which approximate the solutions by means 
of piecewise continuous functions. In this work, we study the DG methods for 
elasticity problem with the derivation of a posteriori error bounds in the L°°(L 2 ) 
norm. Concerning with the study of a posteriori error estimates for elasticity 
system, several different approaches have been developed for stationary prob¬ 
lem: we refer to [6] and [7j for residual type error estimators and method based 
on dual variables, respectively; on the other hand, in [ 7 ] and 13 they consider 
a posteriori error estimates for isotropic problem. In case of time-dependent 
problem, to our knowledge there is still no work for a posteriori error control 
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for elasticity equation, so this is the first study about this. Inspired from the 
recent work of [5] on a posteriori error estimates for finite element method for 
wave equation, we here propose an a posteriori bound in the L°°(L 2 ) norm of 
the error for the elasticity equations, and in case of DG method. The theory 
is developed for both space-discrete case, as well as for the case of fully dis¬ 
crete scheme with the implicit backward Euler method in time discretization. 
The main techniques for this derivation conclude the special testing procedure 
introduced by |9J, and adapt the technique about elliptic reconstruction used 
in [TO] . HU to stationary elasticity (SE) reconstruction. The SE reconstruction 
technique allows to estimate the error of the time-dependent problem via an 
auxiliary SE equation. 


1. Notations and approximation properties 


1.1. Some functional spaces. Throughout this paper, the usual Sobolev norm 
of H m on E C ! d is denoted by || • ||o,e = |MIl 2 (e) , the standard space, 
norm, and inner product notation are adopted, their definitions can be found 
for example in m- 

For the spaces of vector v = (uj)i<j<d, w = (iCj)i<i<d E X{E) d and tensor 
fields <r = € X(E) dxd , we define the following 

scalar products: 

d d 

v w = er : r = ^ cr^r^; 

i =1 i,j =1 


and the norms 

|u| = y/v ■ v, |er| = \Jcr : a. 


In case v is a vector-valued or matrix-valued function, the corresponding term 
||u|| Si £: is defined in a similar manner. Here we use the same notations for the 
norms of vector-valued and matrix-valued. 


Let X is a Banach space, we will make use of the Bochner space L p { 0, T ; X ), 
1 < p < oo endowed with the norm 


u \\lp(0,T,X) 


(fo lMlx) 1/p > 1 < p < oo; 

ess supt e[0)T] |H|A-, p — oo. 


1.2. Triangulation. Let Eh be a subdivision of H into disjoint open sets { K } 
such that = U/fg 7 ^/\, and assume that the mesh is regular , e.g. it has no 
hanging nodes. Denote by hx ■= diam(/l) the diameter of element K and by 
h := mtixxGTh hie the mesh size. 

We assume that the triangulation to be shape-regular , i.e. there exists a 
constant d > 0 such that d := sup^- G -^ ^ < oo, here hx and px denote the 
diameter of K and the diameter of the largest ball inscribed in K, respectively. 

The above assumptions also imply that the mesh is locally quasi-uniform, 
that is, there is a constant k > 0 such that k < , whenever K and I\' 

share a common edge. 
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1.3. Finite element spaces. Given a partition or mesh Th of fi, we define the 
following so-called broken Sobolev space for s > 0: 

H s (T h ) d = {ve L 2 (n) d : v\ K € H s {K) d , K € %}, (1.1) 

and endow it with the norm 

iMU.Tft = ( E hL) 1/2 - 

KeT h 

Consider the DG finite element space 

V h := {v € L 2 (El) d : v\ K € F r (K) d , K € %}, (1.2) 

here P r (lL) denotes the set of polynomials of total degree at most r on K: 

1.4. Averages and jumps. Let £], be a union of all faces (edges) of the trian¬ 
gulation Th, and denote by Eh = E\ U Ej .f with El the union of all interior faces 
(edges) of the triangulation Th, and E^ being the set of all boundary faces. Here 
we will generically refer to any element of Eh as a face, in both two and three 
dimensions. 

Let K + and I\~ be two adjacent elements of Th- Let x be an arbitrary 
point on the common side e = dK + n dK~. For a vector-valued function 
q € H%T h ) d , s > \, let us denote by the trace of q on e from the interior 
of K±. Then we define average and jump at x € e as follows: 

M = \(i* + V); M = (« + -<T). 

If e is a boundary side (e € ), these definitions are modified to 

{M = <?; M = q 

Introducing the function h defined on e € Eh by 

, | _ ( inin{h,K, hx'}, e € E^, e = dK n dK', 
h ' e ~\h K , ee£%, e = dKndn. 

1.5. Auxiliary inequalities. 

Lemma 1.1. (Inverse inequality, cf. [15]) - 1 e t K be an element of Th- Then, 
there exists constant Ci nv independent of h-x, v but depend on the polynomial 
degree r such that 

i) \/v € P r (/\),Ve € dK, ||u|| 0 , e < Ci nv hf^ /2 \\v\\ 0 ,K , 

ii) W € P r (K),Ve € dK, ||Vu • v e \\ 0 , e < C illv / i " 1/2 ||Vu||o,^. 

By |-I j,K we denote the semi-norm of H 3 (K). Then we also have: 

Lemma 1.2. (cf. [15]) ■ There exists a constant C depending only on the shape- 
regularity of the mesh, the approximation order r, and the dimension d 

\ v \j,K < Chtffff 3 \v\i t x, Vu € P r (K), 0 < i < j < 2. 
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Lemma 1.3. (cf. |15| . Theorem 2.6 ). Let K be an element of the triangulation. 
Let v E H S (I\) for s > 1 and let r > 0 be an integer. There exists a constant C 
independent of v and hx and a function v E ¥ r (K) such that 

VO <q<s, ||u - v\\g tK < Ch^ q \v\ s , K ] 

where p = min(r + 1, s ) and e E dK. 

We also have the following approximation properties (cf. jl6l|): 

Lemma 1.4. Let K be an element of the triangulation with hx its diameter. 


Let 
of v 

v E H S (K), let r > 0 be an integer. There exists a 
and hx and a function v E ¥ r (I\) such that 

constant C independent 


VO < q < s,||u - v\\ q! x < Ch^ q \\v\\ S)K , 

s > 0, 

(1.3a) 


VO < q < s,||u - u||o, e < Ch^ 1/2 \\v\\ SjK , 

s > 1, 

(1.3b) 


VO < q < s,\\v - v |i i6 < Ch^ 3/2 \\v\\ s , K , 

s > 3/2, 

(1.3c) 


where p = min(r + 1, s) and e E dK. 


We also have the trace inequality: 

Lemma 1.5. (Trace inequality, cf. (IT]). Let K be an element of the triangula¬ 
tion, and denote by hx the diameter of K. Let e E dK, the following inequality 
holds 

IMIo,aR — C'tr(h K 1 \\v\\l j x + ^aH|Vu||o )A -), Vu E H l (I\), 
here Ct r is independent of hx and v. 

Remark 1.1. Note that analogous bounds as in the lemmas above can be also 
obtained for vector-valued function. 

2. Mathematical model and DG formulation 

2.1. Mathematical model. We consider the stress tensor cr(u) = (iTjj(ri))i<jj<rf, 
which satisfies the constitutive relationship: <r(u ) = C s(u), or equivalently 

d 

VI <i,j < ^ C ijk i£ki{u ), 

k,l =1 

where C = ( Cijki)ijki is a fourth order tensor, independent of t and satisfying 
some symmetry properties: 

Cijkl — Cjikl — Cijik — Cklij- ( 2 . 1 ) 

We assume that the stiffness tensor C is positive definite and piecewise constant 
in fl, then there exists a constant c* > 0 such that 

V(t ij)ij ~i~ 0) 'y ] Cijkakllij ^ c * E V>°- P- 2 ) 

Remark 2.1. With piecewise constant tensor C, there exists a positive constant 
C* such that HCtHo,^ < C** ||' 7 ”||o,iv for any element K E Th, and matrix r E 

L 2 (n) dxd . 
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We consider the equations of linear elasticity problem of finding the displace¬ 
ment vector u = ( Ui(x,t))f =1 at a point x in and at a time t £ [0, T] such 
that 


pdlui 

u = 0 


d d<jjj(u) 
^ j=1 dxj 


fi in (0, T) xQ, i = l,d, 
on dfl x (0, T\, 


(2.3) 


with initial conditions uq £ HQ(Q,) d , and U\ £ L 2 (0,) d : 


j u(-,0) = ito on 0 x {0}, 
( dtu(-,0 ) = u\ on 0 x {0}. 


(2.4) 


Here f(x,t) = (fi(x,t))f =1 € L 2 (0, T; L 2 (0) d ) is a general source function; 
and p(x) is the mass density of the material. For simplicity, we will take the 
coefficient p = 1. 

Assume that the analytical solution u in (12.31) satisfies the following varia¬ 
tional formulation: 


{dt t u, v)n + a(u, v) = l(v), Vv £ Hq (0) <z , a.e. in (0, T), 
where the bilinear form a(-, •) is defined by 


a(u, v ) 


(<r(u),e(v))n 



and the linear form l(-) is given by 


(2.5) 


( 2 . 6 ) 


l(v) 


/ f ■ vdx. 

Jn 


(2.7) 


2.2. DG formulation. 

Definition. Define the bilinear form ah on H s (Th) d X H s (Th ) d , s > 3/2 by 



where v e is an unit normal vector associated to the face e, oriented from K + to 
K~. 

In formulation (12.81) above, the function a in the last terms penalies the 
jumps of u and v over the faces of Th', a = oh -1 with a is a positive parameter 
independent of the local mesh sizes will be specified later to assure the coercivity 
of the bilinear form; a and £ are the stress tensor and symmetric gradient 
respectively, taken element-wise. 

The semi-discrete DG approximation to (ESJ) then reads as follows: Find 
Uh : [0, T] —> Vfr such that 

(■ dt t u h ,v) n + a h (u h ,v ) = (f,v) n for all v £ V h ,t £ (0,T], (2.9) 
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with initial values 


Ufi\t=0 — n h u 0 , 
dt'U'h\t=o — ri/jtti, 


( 2 . 10 ) 


here I - !/, is the orthogonal L 2 -projection onto V } x . and the discrete bilinear form 
ah on \h x V/! is given in (12.8|) . 


Lemma 2.1. (Well-posedness of the bilinear form ah). Let the interior penalty 
parameter be defined as in \2.8 1) . Then there are positive constants k = min{^, 
M = max{2, ‘loT 1 l 2 C- msr C* + 2 C*}, and a m ; n = 4C'? 1V (C'*) 2 c” 1 independent of 
h such that 

\a h (u,v)\ < M\\\u\\\ • |||v|||, Mu,v E V h , 
and for a > a m i n; we have that 

ah(u, u) > k|||u||| 2 , Mu e V/j, 

with DG norm 

lllulll = ( ll e Nllo,A' + “h^^lMllo.e) 172 - 

KeT h e£S h 


3. Semi-discrete error estimates 

3.1. Stationary elasticity reconstruction. In this section, we introduce the 
stationary reconstruction which is the main tool to derive the error of the time- 
dependent problem. 

Definition. (SE reconstruction and error splitting). Let Uh be the DG ap¬ 
proximation given by (12.911 . Let also : L 2 (Lt) d —> V/, be the orthogonal L 2 
-projection operator onto the finite element space V^. We define the SE recon¬ 
struction w = w(t ) € HQ(Tt) d of Uh = Uh(t ) at time t E [0, T] as the solution 
of the SE problem 

a(w,v ) = (g,v)u for all v € #o(D) d , (3.1) 

where 

g := Bu h - 11 h f + /, (3.2) 

and B : V h V h is the discrete operator defined by 

(Bz, v)a = a h (z,v) for all v E V h , 

for each z E V^. 


We decompose the error as 

e := Uh — u = p — 6 , (3.3) 

where 6 := w — Uh and p := w — u E l?d(D) d . 

Remark 3.1. (Well-posedness of w). The SE reconstruction above is well 
defined. Indeed, Buh E is the unique L 2 —Riesz representation of a linear 
functional on the finite-dimensional space V/j. And the existence and uniqueness 
of the solution w in (13.1H . with data Buh — I ~\hf + / E L 2 (fil ) d , follows from 
the Lax- Milgram theorem. 
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Remark 3.2. (The role of w). Consider the SE problem of finding w £ 
satisfying 

- V • (<r(w)) = g, (3.4) 

with g defined by (13.2[) . Let wh £ be the DG approximation to w, defined 
by the finite-dimensional linear system 

a h (w h ,v) = ( Bu h - 11 h f + /, v)n, for all v £ 

this implies a h (w h ,v) = ( Bu h ,v)n = a h (u h ,v ) for all v £ V h , i.e., w h = u h . 
Therefore, the SE reconstruction w is the exact solution to the SE problem (I3.4|i 
whose DG approximate solution is Uh- 

3.2. Semi-discrete error relation. 

Lemma 3.1. With reference to the notation of decomposition as in IIS. 3\) . the 
following error relation holds 

(df t e, v)q + a(p, v) = 0 for all v £ Hq (Ll) d . 

Proof. We have the following expressions 

{dt t e, v)n + a(p, v ) = (8&u h ,v)si + a(w, v ) - (tifcu, v) n - a{u , v) 

= (d? t u h ,v) n + a(w, v ) - (/, «) n 
= (9t t u h , n h v) n + a(w,v) - (/, «)n ^ ^ 

= - a h (u h , V\ h v) + (/, V\ h v)a + a(w,v ) - (/, v) n 
= ~ah{u h , + a(w,v) + (Id h f - f, v)a 
= 0 , 

where in the first equality, we used the decomposition (13.31) , in the second equal¬ 
ity we made use of (12.51) . in the third and fifth equalities, the properties of the 
orthogonal L 2 -projection and the formulation (I2.9|) are used, finally the last 
equality follows the identity ah(uh, — (I— f,v)n = a(w,v), which is 
deduced from the construction of w as in Definition nm □ 

3.3. Abstract semi-discrete error bound. We have the following result for 
controlling the error ||e||^oo(o,T;L 2 (Q) d ) i n terms of the nonconforming error 6: 

Theorem 3.1. (Abstract semi-discrete error bound). Let u and be the 
weak solution in h2. 51) and its DG approximation defined in \2.9j) respectively, 
let w be the SE reconstruction of u^ as in Definition \3.1[ Applying the error 
decomposition {3.31) that p = w — u, 6 := w — Uh, the following error bound 
holds: 


\P\ \L°°(0,T-,L 2 (£l) d ) ^ V^(||Mo-t*h(0)||o,n + ||«(0)||o,n) 


+ 2 


/ ||<9t0||o,ndf + 2 Cfoc* 1//2 ||mi — d t Uh{ 0)||o,n, 

Jo 


(3.6) 


Moreover, we also have 

ll e llL°°(0,T;L 2 (O) d ) ^ II^IIL°°(0,T;L 2 + >/2(||wo-«h(0)||o I n + ||fl(0)||o I n) 

rT 


+ 2 


/ 

Jo 


't^llo.ndt + 2 C'foc* 1//2 ||«i — 3tW/i(0)||o,n, (3.7) 
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where Cpn is the constant of the Poincare’s inequality ||u||o,o < CfoIIV u||o n> 

\/ueH^(n) d . 

Proof. The proof of this estimate is based on the integration by parts and a test 
function which is introduced by [9], has been widely used for example in m, 0 
Firstly, we will prove the bound \3.b\) . Let the test function v : [0,T] X —>• M d 
defined by 

*(*,-)=/ p(s, -)ds, t € [0, T\, (3.8) 

for some fixed t* G [0,T]. Clearly v G Ho{Q) d and p G -ffg(fl) d . Also, we 
observe that 


v(t*,-) = 0, and dtv(t, ■) = — p(t, ■) a.e. in [0,T], 

Set v = v in the error relation m, integrate between 0 and t* with respect 
to the variable t, and integrate by parts the first term on the left-hand side to 
obtain 

-[ (d t e,d t v)ndt + (d t e(t*),v(t*))n - (d t e(0),v(0))n + [ a(p,v)dt = 0. 

Jo Jo 

Due to the symmetry of the tensor C, and using the fact that dtv = — p and 
v(t*, ■) = 0, integrate on [0, T], we rewrite the identity as 

J 0 So ^°( i K i )>*K*)) d * = f o (dtO,p)ndt+(dte(0),v(p))n 

which implies 

^l|p(Ollo,a - ^l|p(°)llo,n + ^a(ti(0),w(0)) = [ {d t 0, p) n dt + (<9 f e(0), «(0)) n . 


/o 


Hence, we deduce 

jimilwi - )||P(0)IIL + io(»(0).i(0)) 


< max 
o <t<T 


||p(i)||o,n / Iloilo,ndf + ||5 t e(0)||o i o||h(0)||o i o. 
Jo 


(3.9) 


From the Poincare’s inequality, the bound of tensor C, and the Korn’s inequality 
(cf. [15|, HU) that ||Vw||o,n < \/2||e(v)||o,n, Mv € we get that 


^ a (f)(0), t>(0)) > ^c*||e:(6(0))||^ a 


1 ^ ri— 2 | 


(3.10) 


> -c*||Vw(0)|| 0in > ^c*C Fn ||n(0)|| 0)a , 
Combining this bound with (13.91) . we arrive at 


llp(Ollo,n - 9llp(°)llo,n < 0 ™ax ||p(t)||o,n / Iloilo,ndt 


(3.11) 


+ ||^te(0)||o,n||w(0)||o,n - —c*C7 Fn ||w(0)|| o>n 
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Moreover, by estimating the last two terms, we arrive at 

^llp(Ollo,fi - ^llp(0)llo,n < max l|p(*)llo,n [ I|3s0||o,nd* + c* _1 C , J. n ||3te(0)||§, 

— — J 0 

(3.12) 

Now, we select t* such that ||p(t*)||o,n = maxo<t<r ||p(t)||o,n- In (13.12|) . using 
the inequality AB — jA 2 < B 2 with A = ||p(i*)||o,f2> B = ||5t^||o,odt and 

then taking square root of the resulting inequality yields 

l|p|lL oo (o,T;L 2 (fi) d ) — V / 2||p(0)||o,n + 2 f ||9tP||o,ndt + 2Cfqc* 1 ^ 2 ||3te(0)||o,Q; 

J o 

(3.13) 

Finally we will get the estimate (13.71) as a consequence of the above estimate. 
Indeed, using the bound ||p(0)|| 0 ,n < ||e(0)|| 0 ,n + ||0(O)|| o ,n; e(0) = u h { 0) - u 0 , 
and dte( 0) = dtUh( 0) — Ui, we conclude that 

ll e llL°°(0,T;L 2 (n) ti ) ^ ll 0 llL~(O,T;L 2 (n)rf) + l|PllL°°(0,T;L 2 (n) ti ) 

^ ll^llL°°(0,T;L 2 (n) d ) + V2(||wo-w fc (O)||o,n + ||e(0)|| Ol n) 

+ 2 f ||9 t P||o,Qdt + 2 Cfqc* 1//2 ||«i - ^^^(0)110,n- 

Jo 

(3.14) 

□ 


Now it remains to estimate the norms involving the conforming error 6 = 
w — Uh by a computable quantity. According to Remark 13.21 we have that that 
lift is exact the DG approximation of w. We have then the following result. 

Lemma 3.2. Given r £ L 2 (VL) d , consider the stationary elasticity problem of 
finding z € ^(D)^ such that 

— V • ( a(z )) = r, z = 0 on 5D, (3.15) 

whose solution can be approximated by Zh £ V h of the following DG method: 

a h (z h ,v) = (r,v)a, forallveV h . (3.16) 

If we assume that an a posteriori estimator Sip exists, i.e. 

\\z - 2 ft||o,n < £ip(zh,r,T h ), 

then we can bound the terms in Theorem as follows 

ll^llL°°(0,T;L 2 (D) d ) < ll<%> ( u h, 9, 7ft) ||L°°(0,T) j (3.17a) 

\/2||0(O)|| o ,n < V2S lP (u h (0),g(0),T h y, (3.17b) 

2 / ||$0|| o ,ndt < 2 [ Sip(d t u h ,d t g,T h )dt- (3.17c) 

Jo Jo 

where g = Buh — Id hf A f , and under the assumption that f is differentiable 
in time. 
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Proof. The first two bounds of (13.1711 are directly derived from the property of 
w as presented in Remark 13.21 For the third estimate, noting that o(-, •) and 
•) are independent of t, then there holds 

a(d t w,v) = (. d t g,v)n , for all v € #o(^) d , 

and 

a h (d t u h ,v) = ( d t (Bu h ),v)n = (d t g,v)n, for all v <E V h , 
noting that here we make use of the fact that the projection 14^ commutes with 
the time differentiation. The estimate is then the difference between dtuh and 
its reconstruction dtw. □ 


4. A POSTERIORI RESIDUAL BOUNDS FOR THE STATIONARY ELASTICITY 

PROBLEM 


In this section, we present two methods to derive a posteriori error bounds 
in L 2 -norm for the stationary elasticity problem. The first one is based on the 
method of duality, and the second one is based on rewriting the DG formulation 
with lifting operator to derive the error estimate in energy norm which will 
provide a bound in the L 2 -norm of the error as well. 


4.1. A posteriori error bound in L 2 -norm by duality method. 


Theorem 4.1. (Error obtained by duality method). Let z € be the 

solution to the stationary problem iS.150 and Zh € be the DG approximation 
of z as in PL 16}) . Then the error for the DG method is estimated by 

\\z ~ Zh ||o,G < &[p{Zh,r,Th), 


with <£jp is given by 

(z h ,r,T h ):=cl Y, h K\\r + V-(<T(z h ))\\l K 

1 K&T h 


+ Y h3 HM^Klllo,e 



where C is a positive constant independent of Zh,r,h and Th.. 


Proof. In this proof we denote by e = 2 - 4 . We consider the adjoint problem 


—V • (a(4>)) = e, in Ll 

<f> = 0, on dLl, 


(4.1) 


and assume that this problem is regular in the sense that the solution 0 satisfies 
0 G i7 2 (fl) rf with continuous dependence on e: 

|| 0 || 2 ,u < C||e||o,o- (4.2) 


This assumption is known to hold, in particular, if the domain is convex 
(cf. [18]). On the other hand, along with Eq. (14.11) . we get the following expres¬ 
sion 
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Now applying integration by parts on each element and using the fact that the 
term [(cr( 0 ))i/ e ]| = 0 a.e. on e, Ve E we obtain 


2 _ 

o,n ~ 


Y / ( CT (0)) : e( e ) d * - X] / 

K&T h JK e&£ h Je 


{{M 0 )K 1 ■ [e]dA 


(4.4) 


Let 0* be a continuous interpolant of 0 as in Lemma fl. 41 satisfying properties 
in m, note that with this interpolation we also have that (f>* = cf) = 0 on 
dQ (see [16j). Thanks to the consistency of the bilinear form ah, we get the 
orthogonality equation a^(e, <j>*) = 0. Then the following expression holds 



noting that the jump terms Ylee£ h f e € ■ [0*Jd A and f e 10*1 ' 

[e]d7l vanish due to the continuity of <j>* which is followed by our choice. 

Subtracting the above equation from (14.41) . by noticing that J K (cr(<f>)) ■ 
e(e)da; = f K (er(e)) : e(4>)dx we have the following equation 



Integrate by parts the first term from (14.61) . we obtain 

Y [ {f+ V ■ ((T(z h ))) ■ ((f>-<j>*)dx- Y l<j>~ 

<eT h jK 

Y [ - 0*)) I/ eJ • {ZhjdA. 


l e llo,n — 


e&£l' 


e&£ h ' 


Then by employing the Cauchy-Schwarz’s inequality, we have that 

l|e||o,n< ( Y h K\\f + V- (< t ( z /i ))||^) 1/2 ( Y h K 4 H ~ <t>*\\o,K) 1/2 

K&T h KeT h 

+ ( Y U3 Wl( a ( Z h)) U ej\\ 2 0,e) 1/2 ( Y h_3 H^ “ ^*llo,e) V2 


e ^ £ h 


e ^ £ h 


(4.7) 


+ ( Y h - 0*)K]f 110,e ) { Y h IMIo,e) 

e&£ h ee£ h 


2 \ 1/2 


Now applying the properties of the interpolation <j>* as in (11.41) . we have the 
following estimates for an element K of the triangulation and an edge or face e 
on OK: 

110 ~ 0*llo,R' ^ Ch%\\ <f>\\ | K 
110 — 0* llo,e ^ Ch 3 K \\m K 

II ( <T (0 — 0*)) l/ e|lo,e < Ch K \\(t>\\l t K 


(4.8) 
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Finally, applying the above approximation results to 114.71) . we deduce 

ll e llo,G ^ hxWf + V • <r(zh)\\o t K 

' K£ 

\ 1/2 (4.9) 

+ h 3 ^||[(^)Klll§, B +E MIMII/J Mk«, 

eg£^ e &£h 

noting that in the last two inequalities, we make use of the locally quasi-uniform 
of the triangulation. 

Then making use of the regularity assumption (14.21) of the dual problem which 
allows to bound the norm of <f> by the norm of e, the proof of the theorem is 
completed. □ 

4.2. A posteriori error bound in i 2 -norm via energy norm. 

4.2.1. Extension of the bilinear form, ah on larger space. To state our a posteriori 
error bounds, we define the space V (h) = 77q (0) rf + V^. On V(/i), we define 
the mesh-dependent energy norm for v G V(/i): 

INI = ( ll e ( v )ll<U + a IINIIo, e ) 1/2 - (4.io) 

K&Th ee£ h 

We shall extend the DG form ah in (12.81) to the space V (h) x V (h) by using the 
lifting operator. This way of expending the bilinear form can be seen from M- 


4.2.2. Lifting operator. We denote by the following space: 

£/» = {T G L\n) dxd : t\k e ¥ r (K) dxd , K G T h }. 

Observing that s(Vh) C and for the piecewise constant tensor C, we also 
have Ce(V^) C 

Definition. For v G V(/i), the lifting operator «£? : V (h) —>• is defined by 

f Jf(v) : rdx = ^ f {vj • |(C r)u e } d A, Vr G 'E h . 

e&£ h Je 

Note that «£? is well-defined, according to the Riesz representation theorem 
for the L 2 scalar product in (see e.g. m, 0 ). 

Lemma 4.1. (Stability of lifting operator). There exists a constant C > 0 
independent of the mesh size such that 

im»)llo,n < C,„ V C*( £ h-'IIHIIg .,) 1 ' 2 

ee£h 


, for any v G V(/i). 
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Proof. By the definition of the operator Jz? and from Cauchy-Schwarz’s inequal¬ 
ity we have for any v € V (/t): 

ll-SfMIkn - sup f a *W-*** 
z ev h Iloilo,n 


= sup 

zeVft 


< sup 
zeV h 


'Eee£ h feM't( Cz ) V ‘& dA 


(2 


||2||o,a 


1/2 


(E eG£h h||{{(C 2 )i/ e i 


2 

0,e 


1/2 


Fllo.n 


< C mv C*( Y, h 

e^Sh 


—1II IT-.H l|2 \ 1/2 




,e) 


O.e) i 


(4.11) 


noting that in the above proof, we have used the inequalities 

h i/2 || {{{Cz)u e }} ||o,e < ^C m v (||Cz||^ + + ||Cz||^_) ' , 

coming directly from the inverse inequality m and the locally quasi-uniform 
of the mesh. □ 


With these definitions, we can now introduce the following (non-consistent) 
DG form on V(/t) x V(h): 


A(u,v) = Y"' j <j(u) : e(v)dx — Y"' f s(u) : Jf(v)dx 


KeT h 


K&T h ' 


Y [ e(v): -S?( u)dx + Y [ a N ' H d -4; 

K&T h Jk ee£ h Je 


(4.12) 


where a is a positive constant which we will consider later on to establish the 
coercivity of the bilinear form A(-, •). 


The bilinear form A above is continuous and coercive on the entire space 
V(h) x V(h). Setting a min = 4Cj 2 nv (C*) 2 c“ 1 . There are constants M and k 
independent of h such that 

|z4(i(, v)| < M|||w||| • |||i;|||, Vu, v € V(/i), 

and for a > a m ; n : 

A(u, u ) > ac|||ix||| 2 , Mu € V(/r), 
with DG norm |||-||| defined in (14.101) . 

Furthermore, since Al = a/ l onV/ l xV/ l , and A = a on Hl{Pi) d x H ( \ (D) d , the 
form A can be viewed as an extension of the two forms ah and a to the space 
V(h) x V(/i). 

The discrete problem then can be equivalently stated as follows: 

Lemma 4.2. (Non-consistent DG formulation). The DG formulation (|3.16l) is 
equivalent to finding Z/, € V/j such that 

A(z h ,v) = (r,v)n for all v € V/j. 


(4.13) 
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4.2.3. A posteriori error bound. The proof of the a posteriori error bound is 
based on the approximation results in US] which allow us to find a conforming 
finite element function which is close to any discontinuous one. 

Lemma 4.3. ("Bounding the nonconforming part via jumps, cf. |19| ). Let V/, 
be the space of piecewise polynomial functions as defined in m, then for any 
function Zh € there exists a function zf £ V c = V/j n HQ(Lt) d such that 

ii*fc- z h\\l,a < eg e niiwiih). 


E ll v (** - z *)ll Ik < C 2 ( E h _1 ll Mile). 

K&T h e&£ h 

where C\.C -2 > 0 constants depending on the shape-regularity of the triangula¬ 
tion. 


Theorem 4.2. (L 2 — a posteriori error bound via energy norm). Let z £ HQ(Q) d 
be the solution to the stationary elasticity problem PL 1 3\) and Zh £ V h be the 
DG approximation of z as in \f.lb | j. Then 

I \z ~ z h \\o,n < &ip{z h ,r,Th), 

where 


&ip(z h ,r,Th) ■= C 


h K\\ r + V ■ (r( z h)\\l K + h III( <T ( Z ^)) I/ elllo, e 

K&T h e &£{ 


+ 


^(h + ah- 1 )||[z,]|| 2 ie 
e££ h 


1/2 


where C is a positive constant independent of Zh,r,h and Th- 


Proof. We decompose the error as follows 

z-z h = z-zf +zf-z h , 

e c €Hg(f2) d e d£Vh 

noting that e c £ i7g(n) rf and e d £ V/,. where z^ £ V c is the conforming 
approximation of Zy, from Lemma 14.31 

Let e c denote an approximation of e c in the space of element-wise constant 
vector functions as in Lemma 11.31 we have that e c £ V/,.. Define r/ = e c — e c , 
we then have that 


A(e, e c ) = A(z, e c ) - A(z h , e c ) = (r, e c )o - A(z h , rj ) - A(z h , e c ). 

= (r,ri)n - A(z h ,ri). 

Moreover, noting that A(zf l ,e c ) = (r,e c )n, then we have 

c *lk(e c )||o,Q < A(e c ,e c ) = (r,?y)n - A(z h ,r)) - A(e d ,e c ). (4.14) 

Now our purpose is to bound the terms on the right hand side of \4-Uy) ■ 
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i) Making use of integrating by parts element-wise A(zh,rj), then the sum of 
the first two terms on the right hand side of (I4.14|) can be written as follows: 


(r,r]) n - A(z h ,rf) = Y [ ( r + V ' a ( z h)) • V^x - Y f M • l(a(z h ))u e }dA 
K&T h Jk e&El Je 

+ Y j K ^ ,<yZh " >: e ^ dx ~ Y ^ dA 


K&T h 


eG £ h ' 


(4.15) 


Now we work on the parts on the right hand-side of (|4.15p . 

For the first part of fl 4-15\ ), by employing the discrete Cauchy-Schwarz’s 
inequality we have 


Y / (r + V • er(zh)) ■ rjdx\ 

'fi-'T. JK 


K&Th 


|2 N 1/2 


< ( Y h2 K\\r + V-cr(z h )\\l K ) 7 ^( Y V 2 Nllo,A') 

KeT h KeT h 

<C(Y ^l|r + V-rr(z /l )||^) 1/2 ||Ve c || 0 ,n ! 

KeT h 

noting the last inequality comes from the approximation (1 1.3 D . 
For the second part of ( | 4.15 ) 

Y [ ivj i(^(z h ))v e }dA\ 

e 


(4.16) 


<^( XI h 1 Mil,8 k) 1/2 ( XI h IIIM z /i)KJIIo, e ) 


2 \l/2 


(4.17) 


K&T h 




Now to estimate the first factor of f!4.17[) . we will use the trace theorem (11.51) 
in conjunction with the approximation property (11.31) : 


h 1 \\v\\o , 8 K < Ctrh 1 {h K 1 \\ri\\l )K + h K \\Vr)\\o,K) 

< C {h K 2 \\r]\\l K 4- ||Ve c ||o iA -) 

< C\\\/e c \\l K ] 


(4.18) 


noting that in the last two inequalities, we make use of the fact that Vry| = 
Ve c \k for all K £ Th, the locally quasi-uniform property of the mesh and 
the inverse inequality (11.21) . then (14.171) infers that 


E 


e ^ E h 


\{a{z h ))u e \AA\ < C'||Ve c || 0 ,A: x ( Y h III( cr ( z /*)) I/ elllo,e) 


2 \ 1/2 


e ^h 


< 


C||Ve c ||o, n ( XI hllKff^KlIlSJ 


2 \l/2 




(4.19) 
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For the third part of \f.l5 ) 

Y jj ^( Zh ) ■■e('n)dx\ < ( Y \\^( z h)\\o,K) 1/2 { Y Mv)\\Ik ) 1/2 


KGT h 


KGT h K&T h 

< h_ 1 HNlllo,e) 1/2 ( Y H e fa)H<U) 

e&£h K&Th 

2 \ 1/2 


1/2 


< c( h-'iiMiiy 

ee£h 


£(e c i o,n; 


(4.20) 


noting that e{tj)\k = e(e c )| k for all K £ Th- 
And the final part of fl 4-15\ ): 

E /i)M - < c(E h-'iiHiiy 1 ' 2 x <e h- i iiizjii5, e ) i/2 

eG £h e ee£h e&£ h 

<C||Ve c ||„, n x (E hyiMiyy 

c€:£h 

(4.21) 


noting that here we use the same approach as in (14.171) . 
ii) It remains now to estimate the term —A(e d , e c ) from (14.141) . We have that 

-A(e d ,e c ) < C* ( E ll £ ( e ^)llo,^)ll £ ( e c)llo,n + ||-S?(ed)||o,nl|£(e c )||o 

K£T h 

/ \ 1/2 (4.22) 

< C|| E (e c )||„, n ( E 11-Mil^ + E h GllMllo,= ■ 

\K£Th e€£h / 

Observing that ||[[e rf ]|| e = |||| e , Ve G S h . We also have ||e(e rf )|| 0 , k < 
C'llVe^Ho,^ for K E Th, and remark that from Lemma 14.31 the noncon¬ 
forming term can be bounded in terms of the jumps of discrete solution. 
We then obtain the following result 

- A(e d ,e c ) < C'||e(e c )|| 0 .o(E h_1 H Ml§,e) 1/2 , (4-23) 

ee£h 

Then, summing up all the bounds (14.161) . (14.171) . (14.201) . (14.211) and (14.231) . 
using the equivalent norms ||Ve c ||on and || £ ( e c)||o,n in ^o(^) rf 5 factorize the 
term Ve c and hnally deviding both sides of the resulting inequality by Ve c , we 
arrive at the energy bound for the conforming part of the error as follows: 


|Ve c || 0 ,n<C4 Y h K\\ r + v ’ <r(zh)\\o,K + Y h WK cr ( z h)) u ' 
K K&T h 


ej| llO.e 




+ Y ah 

ee£ h 


IL 


1/2 


(4.24) 
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Now to derive the error in L 2 (fl)-norm, we will use triangular inequality and 
Poincare’s inequality. By triangular inequality we have that 

\\ z — z h \\o,n < ll e c||o,n + ||e<i||o,n. (4.25) 

Here, by using Poincare’s inequality for the function e c € we have that 

||e c ||o,n < C'FnllVedlo^. 

Then making use of the bound for the energy norm (I4.24p . we derive the bound 
for L 2 -norm of e c , and combining with Lemma [4.31 to bound the term ||ed||o n, 
we complete the proof. □ 

Remark 4.1. (Residual-based a posteriori error indicators). Notice that the 
sum of the above quantities can be decomposed into element-wise contributions. 
For example in case of the global estimator in Theorem 14.21 we can use the local 
error indicator: 

V k = 

\Jh 2 K ||r + V • (cr(^))||§ ^ + h||I(cr(^))i/ e l||2^ van + (h + ah- 1 )!! \z h \ || 2 ^, 

and note that these contributions provide local measures of the magnitude of the 
residual and are able to provide refinement indicators for adaptive strategies. 
And in practice, because the mesh is taken with small hx, then the term h is 
very small in comparison with ah -1 , so this term (h + ah^ 1 )|||z/ l ]||Q g K can be 
considered as ah -1 1 | \zh\ ||q 9k . 

Remark 4.2. We present here two alternative ways to derive a posteriori error 
bunds for stationary elasticity problem. In te first approach, we make use of 
the duality technique, which is well-known for deriving the L 2 —norm of the 
errrors in finite elements methods. But notice that by this method, we need 
the regularity assumption of the adjoint equation. In the second approach, we 
obtain automatically a bound in L 2 — norm from the bound of energy norm, but 
we also recognize that this bound loses an order of the mesh size in comparison 
with the bound obtained by the duality method. 

4.3. Semi-discrete a posterior* error estimates. Now in view of Lemma l3.2l 
and the error bound in stationary case in Theorem 14.11 or 14.21 we obtain the 
following a posteriori error bound for semi-discrete problem. 

Theorem 4.3. With the notations defined as in Definition \S.l[ assume that f 
is differentiable with respect to time. Then the following error bound holds 

11 ® 11 L°° (0,T,L 2 (f2) d — H < ^Ip( u /iifl , )7h)llL oo (0,T) 

+ 2 / Sip(dtUh,dtg,Th)dt + V2£ip(uh(0),g(0),Th) (4.26) 

J 0 

+ V2K — w/j.(0)||o,q + 2Cfqc* 1/2 ||mi — dtUhfif ) ||o,r 2 ■ 
with g = Buh — Id } L f + f , and <ojp is given by {HP or & 

5. Fully discrete error estimates 


5.1. Fullly discrete formulation. 





18 


LUONG THI HONG CAM AND DAVEAU CHRISTIAN 


5.1.1. Space-time meshes and spaces. We now discretize the time derivative by 
using the backward Euler discretization with a subdivison of the time interval 
(0,T) into subintervals (t n ~ 1 ,t n ),7i = 1,... ,1V, with t ° = 0, and t N = T, and 
with the uniform time-step t = t n — t n ~ 1 . Associated with the time subdivision, 
let T}?, n = 0,..., N be a sequence of meshes which may be different from Tu 1 
when n > 1. We also assume that the sequences of meshes TS , n = 0,... ,1V 
are regular and shape regular, uniformly in n. 

For each n = 0,..., N, we denote by VT a DG finite element space of fixed 
degree r built on the partition 7)” 

V£ := {v € L 2 (Q) d : v\ K € P r {K) d , K € 7^}, 

and denote by £J) the union of all edges (or faces) of the triangulation TS. 


5.2. Fully discretization formulation. The fully discretization reads as fol¬ 
lows: For each n = 1,..., N, find u ^ € VJJ such that 

(d 2 ul,v)n + a%(ul,v) = ( f n ,v)n for all v € V£, (5.1) 

where the bilinear form aJJ is given by 


®h(w.v) = X] . a ( u ) : £ ( v ) dx -^2 {{( ct NKJ • HdA 

KeT n j K e&g n Je 

- • H d ^ + [ a H • M d A 

Je ee£j; Je 


(5.2) 


and f n := f(t n , •), the backward finite difference scheme for the discrete tem¬ 
poral derivatives 


d 2 u n h : = 


du% — du 


n —1 
h 


(5.3) 


with 

dul := for n = 1,2,..., N, (5.4) 

T 

and the initial values (n = 0 ) are given by: 


u h := n£« 0 , du° h = ngui, 

here is the orthogonal L 2 -projection onto the finite element space (al¬ 
though other projections onto can also be used for setting the values of vP h 
and du° h ). 

It’s known that the backward Euler scheme is unconditional stable, from the 
coercivity of the bilinear form a)J(-, •), we have that the problem (15. ID admits a 
unique solution (uJJ)i< n <jv at each time step. 


5.3. Space-time reconstruction. With the notations of finite element meshes, 
spaces introduced in Section 15.21 for fully discrete formulation, we have the fol¬ 
lowing definitions: 

Definition. ( Space-time reconstruction). 
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i) Let , n = 0,..., N, be the fully discrete solution computed by the method m, 
n^; : L 2 (Q) d —> VJJ be the orthogonal L 2 -projection, and B n : 
to be the discrete operator defined by 

for z £ \%,(B n z,v)n = a%{z,v), Vv £ V£. (5.5) 

We define the SE reconstruction w n £ HQ(p,) d of u 7 ^ to be the solution of 
the SE problem 


a(w n ,v) = ( g n ,v) n , Vu € #o(D) d , (5.6) 

with = B n Ufr — l~l T £f n + f ; where l"IJJ be the orthogonal L 2 - projection 

onto the space V£, f° = /(0, •), /” = t” 1 J)t_! /(t, -)dt, for n = 1,..., JV. 

ii) Define the reconstruction dw° £ fi) d to be the solution to the SE 
problem 

a(dw°,v) = (dg 0 ,v)n, \/v £ #o(D) d , (5.7) 

with 

d< 7 ° :=S°(^)-n 0/ 0 + /°, 

recall that 11 denotes the orthogonal L 2 projection onto the space and 
du° h = n°U! as introduced in Section 15.21 

iii) The continuous time-extension ujv : [0, x D -A of is defined 

by 


u N (t) : = 


t — t 


.71—1 


-< + 


f n -f 




(t-t n - i )(t n -t)^ 


■9 2 u n 




(5.8) 


for t £ (f n_1 , t n ], n = 1,..., N, with d 2 u'f l given by backward Euler scheme (15.31) . 
We note that u jv is a C 1 function in the time variable, with UN(t n ) = u 1 ^ 
and dtUN{t n ) = for n = 0,1,..., N. 

iv) The time-continuous extension w jv : [0, X !1 A l (i of is given 

by 


w N (t) : = 


t - t 


n— 1 


-iiL 1 + 


f n -t 


-in 


n—1 


(t-t n ~ l )(t n 

T 



(5.9) 


for t £ (t n ~ 1 ,t n ],n = 1, with d 2 w 11 given by the backward Euler 

scheme as in (15.3p . noting that dw° is defined as in (15.71) . By construction, 
this is also a C 1 function in the time variable. 


Definition. (Splitting the error). Decomposing the error as follows: 


ejy := mjv — u 
'■ = PN ~ @N, 

where 6jy := — u n and p N := — u £ L7g(D) rf . 


(5.10) 
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5.4. Fully discrete error relation. 

Theorem 5.1. (Fully discrete error relation). Under the notation as in Defi- 
nition HO1 for t € t n ], n = 1,N we have 

(dt t e N ,v ) n + a(p N , v ) =((l - n %)dfau N , v) a + p n {t)(d 2 u n h , n)»n 

+ a(uiAr - w n ,«) + (f n ~ f,v) o, /or all v € #o(fl) d ; 

(5.11) 

where I is t/ie identity mapping in I?{Ul) d , and p n (t) := —6r(t — * n +*" 1 ). 

Proof. Noting that d 2 t U]y(t ) = (1 + p n (t))d 2 u^ for t € (t n_1 , i n ],n = 1,...,IV; 
then for all a G iLg(fl) rf we have: 

(<9 2 t ejv, v) n + a(p N , v ) = (<9 t 2 mat, u) n + a(tipv, v) - (/, u)^ 

= ((I - njj)$t«tf, w)n + (^mat, n» a 

+ a(w N ,v) - {f,v)n 

= ((I - n h)dt t u N , v) n + p n (t)(d 2 ul H n h v)n 
~ n lv) + a(w N , v) + (n lf n - f, v)n 

= (0 - n h) d tt u N, «)n + p n (t)(d 2 u%, n h v h 
+ a(w N - w n , v) + (. f n - f, v)q] 

noting that the first equation follows from the formulation (12.5p of the weak 
solution u. In the second and third equalities, the properties of orthogonal L 2 - 
projection on space VJJ are employed, and we also use the formulation of 
fully discrete scheme in the third equality, and the last equality comes from the 
identity a n h {u r f i , njjv) - (n £/ n - v) n = (B n uf. - (n \f n - v) n = 

a(w n ,v), \/v € i7g(fl) d . □ 

Remark 5.1. ( Zero-mean value property). The particular form of the remainder 
p n (t) satisfies the zero-mean value property 



which is an important point for estimating the a posteriori bound derived 
from the term p n {t){d 2 ^,^ v)q in the error relation m above (see later 
in Lemma 19.41 for the proof of time-reconstruction error bound). 

5.5. Abstract fully a posteriori error bound. To analyze the error, we 
introduce the following quantities that will provide error estimator of the fully 
discrete scheme in Theorem 15.21 These indicators are originally introduced by 
Lakkis et al. in |5J concerning a posteriori error estimates for finite element 
method to the wave equation. 

Definition. (A posteriori error indicators). We define 
i) The mesh change indicator is given by 

N N—l 

Cue :=E / 11(1 - n h) d tUN\\o,ndt+^2(t N -t n )\\(r\l +1 - n£)0ujj||o,n- 

n — 1 J I™ 1 n — 1 
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ii) The evolution error indicator reads 


fi N 

Cevo := / Iloilo,ndt, 

J 0 


where Q : (0, t N ] —>• with Q\^ t „-i t „] := <5 n , n = 1,..., N and 


G n (t) := * )2 dg n - 




4r 


W-7n, (5-12) 


with £/ n as in Definition 15.31 and 7 n := 7 n _i + (r 2 /2)dg n + (r 3 /12)d 2 g n , 
n = 1,..., N with 70 = 0; 

iii) The data error indicators are given by 

-\ N / P n _ \ V 2 


which can be viewed as an error estimator related to the time-oscillation of 
the source term. 

iv) The time reconstruction error indicators 


—me 


! n^ 


Ui 




1/2 


Theorem 5.2. (Abstract fully discrete error bound). Let u be the weak solution 
to 112.5 1) : uj\r and are reconstructed from the fully discrete solution { u]t}n=i 
and its SE reconstruction {w n }^ =1 , respectively, as in m and m, and re¬ 
calling the notation of indicators in Definition ] 5. 51 then the following a posteriori 
error estimate holds 

fi N 

ll e Ar||L°°(0,t JV ;L 2 (n)' i ) < \\^n\ \L°°(0,t N -,L 2 (£l) d ) + V^H 0 iv(O)llo,n + 2 / II dtd iv||o,fidi 

Jo 

+ 2 (Cmc + Cevo + Cose + Cr.Rec) 

+ V^H^o — w°||o,si + 2(7 fqc* 1 ^ 2 ||wi — 9w°||o,a, 

(5.13) 


where Cpn is the constant of the Poincare inequality, and du ^ are the or¬ 
thogonal L 2 —projections onto the space defined in Section \5.2i 

Proof. The derivation of this error bound is similar to in senri-discrete case. 
Firstly, we will take a test function v = vjy with ujv is defined similarly as 
in (13.81) 

v N (t,-) = J p N (s,-)ds, t e[0,t N ], (5.14) 

assuming that f m_1 < t* < t m for some integer m with 1 < m < N, and Pat is 
defined as in (15.3p . Then vj\r E HQ(Q) d and p N E ^(D) 6 *, and we observe that 

v N {t*,-) = 0, and d t i>N(t, ■) = -p N (t, •) a.e. in [ 0 , t N ], 
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We integrate the resulting equation with respect to t between 0 and t*, to arrive 
at 


(dt t e N ,v N )nd t+ a(p N ,v N )dt = !>(*•)> < 5 - 15 > 


i —1 


where 


m— 1 r £n 


m*) ■= E / (0 - n Z)$ t u N ,v N hdt+ / ((I - U^)dlu N ,v N )ndt 

it"- 1 it™- 1 

m-l .fn. t * 

l 2 {t*) :== E / cl(wn — w n ,v jv)dt + / a(w]y — w m ,visr)dt 
it' 1 - 1 it"*- 1 

m-i - t n ~ ~ 

EO := E / (f n - f,v N ) n dt+ (f m - f,v N )ndt, 


n=l 
m-l .*■ 


Pt 

/ ^(5 2 <,n^ 7V ) n dt+ / /i m (9 2 <,n^) c dt 

Z ' ] fn — l I fm — 1 

n =1 171 


(5.16) 


Integrating by parts the first term of the left-hand side of Equation (j5. 15|f . , and 
using the properties of i>]\r, we obtain 

i-t* 

rat . 


rt *id 2 rid 
iPJvllo.n^ — 


to 


2 dt 


io 


2 dt 


i(v N ,v N )dt= / ( d t 0 N ,p N ) n dt 

Jo 


+ (<9t e At(0), Vn(0))h + EE(^*)’ 


Z=1 


(5.17) 


which implies that 


IIPJv(**)llo,n “ 2 IIpai(O)|| o,n + ^(EKO), vjy(O)) = (d t 6 N ,p N )ndt 


+ (d t e N (0),v N (0)) n + 


i= 1 

(5.18) 

On the other hand, with the indicators Cmc> Cevoj Cose and CT.Rec introduced in 
Definition 15.51 we have the following bounds: 

Lemma 5.1. The following inequalities hold: 

Zi(0 < CMG o mmc T ||piv(i)||o,n; 

^2(0 < Cev 0 mmc T ||p A r(t)||o,a; 

-^(^ ) Cose ||PAt(i)||o,c; 

^(O < CT.Rec max IIPAT(i)||o,n; 


(5.19) 
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The proof for the above lemma follows from the applications of integration 
by parts, the commutation of orthogonal L 2 -projection with time differentiation 
and time integration, and the zero averages of p n and f n — f n on the interval 
[t n ~ 1 ,t n ], These results for wave equation were presented in f5j, but for the 
sake of completeness, we also present, even they are similar, the proofs for our 
elasticity problem in Appendix [9] 

Therefore we have 

^llminilw! - i|lwK0)|lS,n + ia(tw(0),two)) 

||dt0jv||o,rcdi + Cmc + Cevo (5.20) 

T Cose + Cr.Rec^ + ||^t®iv(0) 11 o,q || ^iv(O) 110,0 • 

We select t* such that ||pjv(i*)||o,f2 = max 0<t<t jv ||pAi(i)||o,f2- Then following 
analogous arguments as in the proof of the semidiscrete case, we end up to 

\\ e N\\L°°(0,t N ;L 2 (n) d ) < II||L°°(0,T;L 2 (f2) d ) + ^2(ll w 0 — w hl|o,fi + II#y( 0) ||o,o) 

+ 2^^ ||^||o,ndt +CMC + Cevo 

+ Cose + Cr.Rec^ + 2 CfQC* ^Hui — 

(5.21) 

The proof is then completed. □ 

6. Fully discrete a posteriori error estimates 

To arrive at a practical a posteriori bound for the fully discrete scheme from 
the abstract error estimate in Theorem 15.21 above, it remains to bound the 

rt N 

terms ||0 jv(O)|| 0 ,r2, l|0Y||z / °°(o,P v ;L 2 (Q) d ) and J Q \\d t 6 iv||o,odt. We will establish 
the bounds of these terms in Proposition 17.11 17.21 and 17.31 below. This enables 
us to prove the following error estimate: 

Theorem 6.1. (Fully discrete a posteriori bound). With the same hypotheses 
and notations as in Theorem El we have 

W^-N lll / oo (0,t JV ;L 2 (f2) li ) < Csp + Ctp + ClC, (6.1) 

where Csp mainly accounts for the spatial error, Ctp mainly accounts for the 
temporal error and Cic represents the initial conditions of the problem. They 
are given as follows: 

Csp = Csp,l T Csp,2 T Csp,3) 

Ctp = 2 (Cmc + Cevo T Cose + CT.Rec) > 

Cic = V2IK — u° h ||o,q + 2CpfiC* 1/2 ||ui — d«°llo,n> 


< max ||pAr(t)||o,n 
0 <t<T ’ 


( 6 . 2 ) 
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where the temporal indicators CmCj Cevo Cosc> CT.Rec are defined in Definition 15.51 
and the spatial indicators Csp,i; Csp,2; Csp,3 are given by 

Csp,i = V2^>, 

Cs P ,2 = max|^^ip(3w^aflf°,7^)|,3 o max v + 2C 2 Fn c~ 1 \\f n - / n ||o,n), 

N N 

Csp,3 = ^(*^ip + ^tp X ) + XT ^TCppc* 1 ll<9/" — 9f n ||o,n, 

n =1 n =1 

(6.3) 


u'i 


where Sffi := <S\p(u^, B 
defined from ( fpy or 


n";/ n + f n ,Tff), for all 0 < n < N, S IP being 


7. Proof of Theorem urn 

This section is devoted to the proofs of Proposition 17. 11 17.21 and f7. 31 Plugging 
these results into Theorem 15.21 we will get Theorem 16.11 
At first, we clearly have the following bound: 


Proposition 7.1. The following estimate holds 

\/2||0jv(O)||o,n < Cs P ,i) (7.1) 

where £ sp! i is introduced in Theorem fOl 

Secondly, we establish an estimate for the term ||^||X/ oo (o,* JV ' ; z. 2 (r 2 j rf ) as f°ll° ws: 
Proposition 7.2. The following estimate holds 


\\ON\\L°°(0,t N ;L 2 (n) d ) < Csp,2, 
where Csp ,2 is introduced in Theorem \6.1\ 


Proof. According to the construction of u n and w n in Definition 15.31 for t € 
(t n , t n ], n = 1,..., N, we have the following expression 

d N {t) = w N - u N 


t-t 


.71 — 1 


t n ~t< n—1 ,.n-U 1 )( tn ~ t ) ,o2...n fl 2„.n 


(w n - ul) + —— (w n ~ v - u n ~ l ) - 


(^v-ax), 

(7.2) 


which yields 

noting that 


I0iv(*)||o,fi < 3 max \\w n - u n 


0 <n<N 


h II 0,f2 1 


max 

te{t n ~ 1 i t n ] 


(t-t n ~ 1 ){t n -ty _ 4t 2 
r “ ~27 


Now we need to estimate the terms ||iu n — Wt||o,a,0 < n < N and ||— 
c?w^||o,n. To estimate the term ||ru n — w)M|o,q, we will prove the following result 
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Lemma 7.1. 


The following estimate holds for all 0 < n < N 
\\ wTl ~ u h\\o,n < ^tp + 2Cpqc^ 1 II/'' — /" ||o,f 2 ) 


(7.3) 


with (off, := £ip(u 
defined from (B 


n T>n„,n 
h B u h~ 


or 


B- 


FI/)/ 71 + f n ,Tjf), for all 0 < n < N where S’ip is 


Proof. To prove the result in Lemma [7. 11 we proceed as follows: first, we define 
w_ n € to be the solution to the SE problem 

a(w n ,v ) = (B n ul - n ir + f n , v)n (7.4) 


for n = 0,1,..., N. 

Note that since f° = f°, we have w° = w °- On the other hand, we know that 
uf is the DG approximation in V)( of the SE problem (|7.4[) . noting that to prove 
this we follow the same proof as in Remark 13.21 Then in view of Theorem 14.11 
this implies that 


||m n - K ||o,a < c$p(uiB n vZ - n n h r + f n ,Vf), 

for n = 0,..., N. 


(7.5) 


Second, we need to estimate \\w n — ru n ||o,Q- Observing that w n — uf n is the 
solution of the stationary elasticity problem with the load f n — f n : 

a(w n -w n ,v) = (f n -f n ,vh , 

we end up to 

\\w n -w n \\o,u < 2C'| n c/ 1 ||/ n -/ n ||o,a, forn = 1,..., N; (7.6) 

due to the stability of SE problem. Finally, thanks to the triangular inequality 

\\w n - uf || 0 ,Q < || w n - ru n || 0 ,f 2 + \\w n - uf || 0 ,Q, (7.7) 

along with the bounds (17.51) . (17.61) imply Lemma 17.11 □ 

Now it remains to estimate ||<9-u ; 0 — <9tt9||o,n. Similarly, we have 

a°h(9u° h ,v) = (B°(du° h ),v) n = (B°(du° h ) - 14 ° h f° + f°,v)n , for all v € V“. 

Hence comparing with the construction (15.71) . du° h is the DG solution of dw°, 
which yields 

||dw° - Su^|| 0) n < gip{du° h ,dg° ,7^). (7.8) 

The proof of Proposition 17.21 is thus completed. □ 

Proposition 7.3. The following estimate holds: 

f t N 

2 / ||<9 t £/v||o,ndt < Csp,3 5 

Jo 

where Csp ,3 is defined by (m- 

Proof. Similarly to the proof of the previous proposition, the construction of 
ujy and wn in Definition 15.31 and (17.2|) . for t € (t n ~ 1 , t n ], n = 1 ,N, implies 
that 

d t e N = dw n - dui + T~\t n - t)(3t - 2t n_1 - t n )(d 2 w n - d 2 ui ), 
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from which we deduce that 
r t n 


It n 


-i 


9i0Ar||o,odf < r\\dw n - dvlf || 0) n, 


noting that 


rt n 

/ {t n — t)(3t — 2t n ~ l — t n )dt = 0. 
Jt "- 1 


Summing up (17.91) for n = 1 ,,N, we get 



N 

||3i07v||o,odt < ^T||<9m n 

72—1 


du n 


iillo.n- 


(7.9) 


(7.10) 


(7.11) 


It remains to estimate the terms \\dw n — du 1 ^ ||o,n- We will make use of 
the triangular inequality by combining the bounds for ||chn" — dull ||o,G2 and 
\\dw n — 3m n ||o i o. From the definition of {w_ n }n=o i n (17-41) . we obtain that 
dw n G i7g(n) rf , n = 1,..., N is the solution of the SE problem 

a(dw n ,v) = ( d(B n ul) - 0(n ID + df n ,vfo, Vv G (7.12) 


We then have that dw n — dw n is the solution to the SE problem 

a(dw n - dw n , v) = ( df n - df n , v) n , Vv G H^(Q) d , for n = 1,..., N. 


Then from the stability of steady-state elasticity problem, we have 

\\dw n - 3m n || 0 ,n < 2Cp n c~ 1 \\df n - df n \\ 0 ,n, for n = 1,..., N. (7.13) 

Now we derive a bound for \\dw n — du?\\o } n in the following lemma. 

Lemma 7.2. For n = l,...,jV, let dw n G iLo(fl) rf be the solution of prob¬ 
lem \7.1ty) and du^ is defined by the backward Euler scheme (5.f) where u^ is 
the fully discrete solution from. H5.1\) . The following error bound holds 

\\dw n - du n h || 0 ,n < (7.14) 

with £* := S’ip( v*,B n v% - + f n ,Tff), for all 0 < n < N and S IP is 

defined by CP9 or 

Proof. We again denote by the solution in ()7.4|) . for i = 0,1 ,,N. From 
the backward finite difference scheme, we obtain the bounds as follows: 


)w n - du% || 0 ,n < 


1 


\w n - «ft|| 0 ,n + Ww 71 - 1 - <; _1 ||o,n) 


(7.15) 


for n = 1,... ,N. Employing (17.51) allows us to end the proof for Lemma [7.21 □ 

The proof of Proposition 17.31 is completed. □ 


Remark 7.1. In case of stationary mesh (i.e. the same mesh is used between 
the initial time t° and the final t N , or Tfi = for all n = 1,... ,N ), we 

have an alternative result for Lemma O as follows: 

\\dwfi - d<|| 0> n < <%> (dul d(B n u n h ) - d(U n h f n ) + df n ,T h n ), (7.16) 
for all 1 < n < N and £\p is defined from ED or ED. 
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Proof. Indeed, from the definitions of {t« n }^Lo hi (17.41) . we deduce for all n = 
1,..., N that: 

a(dw n ,v) = (d{B n u n h ) - d(n m + df n ,v) a , Vu € (7.17) 

On the other hand, from the definitions of the discrete operators {B n }^ =() 
in (ESD and property of orthogonal L 2 -projection, we have that for all n = 

<{dul,v) = (d(B n K),v)n, Vu e V£ 

= o d{B n u n h ) - d(nm + df n ,v) n , v« € vi '' 1N; 

which implies that du is the DG approximation in V)) of the boundary value 
problem (17.171) . so we obtain the error estimate (17.161) . □ 

8. Conclusion 

In this chapter, we have carried out an a posteriori error analysis for the sym¬ 
metric interior penalty Galerkin method for the fully discretization of the time- 
dependent elasticity equation. The work is inspired by the method in [5j which 
is used for the wave equation in case of conforming finite element method and 
is expanded to our elasticity problem in case of DG method. This method com¬ 
bines the SE reconstruction technique, the special testing procedure introduced 
by 13 , and a suitable space-time reconstruction allows to derive a posteriori er¬ 
ror estimate for the time-dependent problem from the error of the auxiliary SE 
equation. We stress that this strategy can be adapted by many DG methods, as 
long as there exists a posteriori error estimate in L 2 — norm of the corresponding 
DG method for stationary problem. The numerical implementation of the pro¬ 
posed bounds in the context of adaptive algorithm strategy will be considered 
in upcoming works. 


9. Proof of Lemma 15.11 

In this section, we present the proofs of Lemma 15.11 in order to prove The¬ 
orem m This proof is similar to the proof carried out for the wave equation 
in [5]. 

Lemma 9.1. (Mesh change error estimate). Under the assumptions of Theo¬ 
rem. \5.2\ and with the notation in 15.1611 . we have 

< CMC^ax^ ||piv(i)||o,n- 

Proof. Observing that the projections ll)), n = 1,..., N commute with the time 
differentiation, we integrate by parts with respect to t, ending up to 

m ~ l pt 71 ft* 

-=y2 (0 “ ^h)9tu N ,PN)ndt.+ / ((I - n™)d t u N ,p N ) n dt 

Ji™- 1 Jt ™- 1 

m—1 

+ E« n r‘ - ni)dut,v N (f‘))n - ((1° - r©ftuS,* w ( 0 ))„. 

71=1 


(9.1) 
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The first two terms on the right hand side of (19. 1 [) are bounded by 


(m— 1 


x rt rt 1 

V / ||(l - n2)$ujv||o,ndt+ / ||(l - n™)d t u N \\ 0t ndt 

n=l Jt™- 1 j 

(9.2) 


Recalling the definition of and that of dtU]\r(t n ) = du n = 0,1,... , IV, we 
can bound the last two terms in the right hand side of (19.111 by 


/ m—1 


max ||pAr(t)|| 0 ,G 
0 <t<T ’ 


51 (** - ‘“HKnr 1 - n i)S<llo,n +1* 


- nj)8u° 


hllo,n 


^ n=l 


(9.3) 

Noting that (I — V\ a h )du a h = 0 because du G h € V°, and bounding 19.31 with t* = T , 
we obtain the result of Lemma 19.11 □ 


Lemma 9.2. (Evolution error bound). Under the assumptions of Theorem \5TA 
and with the notation < 15.16]) . we have 

Mt*) < Cevo^nax^llpAr^Ho.G. 

Proof. First we observe the identity 

w N - w n = -( t n - t)dw n + (r~\t n - l) 3 - ( t n - t 2 )^d 2 w n (9.4) 

for any t in (t n ~ 1 ,t n ], n = 1,... ,m. Hence from Definition 15.31 we deduce that 

a(w N -w n ,v N ) = (-( t n -t)dg n +( K T-\t n -tf-{t n ~t 2 ))d 2 g n : v N )u . (9.5) 

Then the integral of the first component in the inner product on the right-hand 
side of (19.51) with respect to t between t n ~ l and t n is given by Q n (see (15.121) . 
Hence, integrating by parts on each interval (t n ~ 1 ,l n ] gives rise to 

m . t n t * 

h (0 = E/ (G n ’PNhdt+ (G m ,p N hdt, (9.6) 

Jt n ~ 1 If ™- 1 

which implies the expected result. Note that the choice of the constants "y n in 
(15.121) makes G continuous at t n , n = 1,..., IV, and we also have 1/(0) = 0 . □ 


Lemma 9.3. (Data approximation error estimate,). Under the assumptions of 
Theorem El and with the notation \5.1b\ we have 



X 3 (t*) < Cosc ( m^ 1 ||piv(t)||o,n- 

(9.7) 

Proof. We begin with observing the zero-mean value of f n 

— / on [t n_1 , t n ] as 

follows 

rt n 



/ (r - /)dt = o 

Jt n ~ i 

. ,m — 1. Hence we have 

(9.8) 

for all n = 1,.. 



m—1 pfri m—1 pfn 

E/ (/ n - f,v N )ndt = E / (f n - f,v N -v n N )ndt, 

n=l n=l 


(9.9) 
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where v^(-) := r 1 f t n-i VN(t, ')dt. Using the Friedrich-Poincare’s inequality 
with respect to the variable t we have 


i-f _ r 2 H n 


o,n 


dt 


(9.10) 


and recalling that dtVN = ~PN, we get 


m—1 iri—± / 

n=l Jt n=1 


m—1 , r fr\ 


- f llo,nd* 


1/2 / f t n 


\v N — vjvllo.ndi 


n n ~ 


i m—1 / r t n \ 1 / 2 / r tn 

* E (/„- lir - /l| U di ) - 2 iimiio,» df 


< 


2 7T 


m -1 / rt 

3 / 




n=l 


\ 1/2 


I \\f n ~ f\\o,n dt j ||pjv(t)|| 0 ,n. 

(9.11) 


For the remaining term in Z 3 , we first observe that 


/*t* /»£* «£* 

/ ||^ 7 v||o,n di < / W IIPY|lo,odsdt < r 3 max ||p iV (s)||o,n> ( 9 - 12 ) 

J t m-1 y t m-l Jt 0<S<T 


which implies that 


U* _ //•**„ \ 1/2 

/ (f m ~ f,v N ) n dt< ( T 3 \\f n - f\\o tn dt ) max ||piv(t)||o,n, 

Jtm-l \ Jtm-1 J 

(9.13) 

and gives Lemma 19.31 

□ 


Lemma 9.4. (Time-reconstruction error boundj. Under the assumptions of 
Theorem EJ and with the notation L5.16\) . we have 


T 4 (t*) < Ct-Rbc^wc, ||pj\r(f)||o,n- (9.14) 

Proof. Noting that <9 2 m)( is piecewise constant, and recalling the zero-mean 
value (15.111 of p n in [t re_1 , t n ], n = 1,..., IV, we have 

m—1 m—1 pffn 

£ / p n ( 5 2 <, n^jv)ndt = £ / p n (a 2 <, - ^)) n dt, (9.15) 

n=! n=! 
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where v N (-) := r 1 fj'a-i vy(t, -)dt. Since the projections n = 1 
commute with the time integration, we obtain 


m-l „ t n _ 

Y / ^ n {d 2 ul,Ul{v N -v n N )) n dt 

n= 1 


< 


2vr 


m —1 

£ v,, 

n=1 x 
m— 1 / 


\ !/ 2 / ft 

y n d 2 K\\i n dtj < i - 2 

1/2 


t 2 ||njJpAr|lo,n 


< 


m-i , \ 1/2 

7^~Y \ T 3 ||// n 3 2 u^j|o ) max ||p/v(#)||o,G- 

27T z —' V /in-1 / 0<t<T 

n =1 v r 7 - 

(9.16) 


For the last term in Z 4 , by using an argument similar to (19.1211 . we have 

r (p m a 2 <,n^ JV ) n dt< 

Jt rn ~ i 


and this completes the proof. 


t* \ 1/2 

T 3 ||p m 5 2 u^||^ a dt) max ||p/v(t)|| 0 ,n, 

71—1 / 

(9.17) 


□ 


References 

[1] P. Houston, I. Perugia, D. Schotzau, A posteriori error estimation for discon¬ 
tinuous Galerkin discretizations of H (cM?'Z)-elliptic partial differential equations, 
IMA Journal of Numerical Analysis (2005), pp.11-25. 

[2] F. Karakatsani, C. Makridakis, A posteriori estimates for approximations 
of time-dependent Stokes equations, IMA J Numer Anal (2007) 27 (4): 741-764 
2006. 

[3] E. H. Georgoulis , O. Lakkis , J. M. Virtanen, A posteriori error control for dis¬ 
continuous Galekin methods for parabolic problems, SIAM J. NUMER. ANAL. 
Vol. 49, No. 2, pp. 427- 458, 2011. 

[4] Bernardi, C. and Siili, E., Time and space adaptivity for the second-order wave 
equation. Math. Models Methods Appl. Sci. 2005, 15, 199-225. 

[5] E. H. Georgoulis, O. Lakkis , C. Makridakis, A posteriori Loo(L 2 )-error bounds 
for finite element approximations to the wave equation, IMA Journal of Numer¬ 
ical Analysis (2013) 33, 1245 - 1264. 

[6] R. Verfurth, A review of a posteriori error estimation techniques for elasticity 
problems, Comput. Methods Appl. Mech. Engrg. 176 (1999) 419-440. 

[7] Y. Chen, J. Huang, X. Huang, Y. Xu, On the Local Discontinuous Galerkin 
Method for Linear Elasticity, Mathematical Problems in Engineering Volume 
2010, Article ID 759-547. 

[8] T. P. Wihler, Locking-free adaptive discontinuous Galerkin FEM for linear 
elasticity problem, Mathematics of computation Volume 75, Number 255, July 
2006, Pages 1087-1102. 

[9] Baker G. A., Error estimates for finite element methods for second order hy¬ 
perbolic equations. SIAM J. Numer. Anal. 1976, 13, 564-576. 

[10] C. Makridakis, R. H. Nochetto, Elliptic reconstruction and a posteriori error 
estimates for parabolic problems, SIAM J. Numer. Anal. 41 (2003), no. 4, 1585- 
1594 (electronic). MR MR2034895 (2004k:65157). 

[11] Lakkis O., Makridakis C., Elliptic reconstruction and a posteriori error esti¬ 
mates for fully discrete linear parabolic problems, Math. Comp., 75, 1627 - 
1658. 




A POSTERIORI ERROR ESTIMATES FOR DG METHOD TO ELASTICITY PROBLEM 31 


[12] Adams R. A., Fournier J. J. F., Sobolev Spaces: Pure and Applied Mathematics, 
vol. 140, 2nd edn. Amsterdam: Elsevier/Academic Press, pp. xiv+305. 

[13] 1. Babuska and W.C. Rheinboldt, A posteriori error estimates for the finite 
element method, Internat. J. Numer. Methods Engrg. 12 (1978) 1597 - 1615. 

[14] M. J. Grote, A. Schneebeli, D. Schotzau, Discontinuous Galerkifinite 
element method for the wave equation Siam J. Numer. Anal. Vol. 44, No. 6, pp. 
2408-2431. 

[15] B. Riviere, Discontinuous galerkin methods for solving elliptic and parabolics 
equations: Theory and implementation SIAM. 2008. 

[16] Beatrice Riviere, Mary F. Wheeler, Vivette Girault,A priori error estimates for 
finite element methods based on discontinuous approximation spaces for elliptic 
problems SIAM J. NUMER. ANAL. Vol 39, No. 3, pp. 902-931. 

[17] Ohannes A. Karakashian, Frederic Pascal, Adaptive discontinuous Galerkin ap¬ 
proximations of second-order elliptic problems, European Congress on Compu¬ 
tational Methods in Applied Sciences and Engineering ECCOMAS 2004. 

[18] P. G. Ciarlet, Mathematical Elasticity: Three-dimensional elasticity, Volume 

I, Antony Rowe Ltd, East bourne 2005. 

[19] Ohannes A. Karakashian, Frederic Pascal, A posteriori error estimates for a 
discontinuous Galerkin approximation of second-order elliptic problems, SIAM 

J. NUMER. ANAL. Vol. 41, No. 6, pp. 2374-2399 , 2003. 

[20] D.N. Arnold, F.Brezzi, B. Cockburn, L.D. Marini, Unified analysis of discon¬ 
tinuous Galerkin methods for elliptic problems, SIAM J. Numer. Anal. 2002, 
Vol 39, No 5, pp 1749 - 1779. 

[21] T.H.C. Luong, A posteriori error bounds for elasticity equations- Asymptotic 
expansion for Stokes problem, PhD Thesis, University of Cergy-Pontoise, France 
(10/2014). 

Department of Mathematics, University of Cergy-Pontoise, France 
E-mail address: thi-hong-cam.luong@u-cergy.fr 

Department of Mathematics, LTniversity of Cergy-Pontoise, France 
E-mail address: christian.daveau@u-cergy.fr 



